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The dynamics of temperature fluctuations of a gas of Brownian particles in local 
equilibrium with a nonequilibrium heat bath, are described using an approach con- 
sistent with Boltzmann-Gibbs statistics {BG). We use mesoscopic nonequilibrium 
thermodynamics {MNET) to derive a Fokker-Planck equation for the probability 
distribution in phase space including the local intensive variables fluctuations. We 
contract the description to obtain an effective probability distribution {EPD) from 
which the mass density, van Hove's function and the dynamic structure factor of 
the system are obtained. The main result is to show that in the long time limit the 
EPD exhibits a similar behavior as the superstatistics distribution of nonextensive 
statistical mechanics (NESM), therfore implying that the coarse-graining procedure 
is responsible for the so called nonextensive effects. 
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I. INTRODUCTION 

The existence of fluctuations of intensive parameters, such as the inverse temperature, the chemical 
potential, the pressure or the energy dissipation rate in a turbulent fluid, is a common feature in a 
variety of physical systems in nonequilibrium states j^, j^, [J]. However, the theoretical description 
of these fluctuations is well beyond the range of apphcability of the theory of fluctuations as developed by 
Onsager and Machlup Q| and subsequent generalization by Fox and Uhlenbeck [gI, and others I?]. This 
is essentially due to the fact that these theories are restricted to describe the dynamics of fluctuations 
of extensive state variables only. 

One of the most used methods to deal with intensive variable fluctuations is the superstatistics ap- 
proach of NESM, 8], [qI, 11 1 . This theory claims that a proper description of these fluctuations 
for systems with a sufficiently complex dynamics, demands the introduction of a more general statistics 
— Tsallis statistics (TS) — than the usual Boltzmann-Gibbs statistics (BG). Furthermore, according 

n n 

to Refs. [3j, [Uj, TS is only one of many more general statistics — superstatistics — that can deal with 
the physical situations mentioned above. An important and simple system where these issues have been 
explored is Brownian motion in the presence of inverse temperature fluctuations In this reference 
the fluctuations of the intensive parameter 1/T are not described as a fully time-dependent phenom- 
ena, rather only their statistics is taken into account and postulated without a physical justification, 
by assuming that if /? is a distributed random variable. By averaging over its (static) statistical 
distribution this model generates TS for v for a stationary state. On this basis, it is suggested that 
BG is no longer capable of describing correctly this physical situation and that a new statistics (TS) 
is required. 

The main motivation of the present work is to consider inverse temperature fiuctuations in the same 
model as in Ref. iJll, but describing their dynamics in terms of a stochastic process, and following 

□ n n 

their time evolution towards equilibrium. By using MNET [1^, [13[, PJ, the time evolution of the 
nonequilibrium state is described in the phase space of the particles, including the local intensive 
fluctuations as state variables, by means of a Fokker-Planck equation [FPE). In the diffusion limit this 
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description leads in a natural way to a stochastic Smoluchowski equation for the contracted distribution 
unction. By following a method introduced by van Kampen to solve multiplicative stochastic equations 
isl . we derive an equation for the effective distribution function (EDF), which still contains all the 
effects produced by the induced (intensive) fluctuations. The long time limit behavior of the EDF 
as a function of position, is similar to the behavior found in NESM for the so-called superstatistical 
distribution. In this way, our approach shows to be an alternative way to describe the dynamics of 
intensive fluctuations consistently with BG. 

lol. In Section 3 we 



We will proceed as follows. In Section 2, we present a derivation of the FPE, 
model the temperature fluctuations by means of a stochastic process and after contracting the FPE, 
we derive a stochastic multiplicative Smoluchowski equation (MSE), valid for an arbitrary stochastic 
dynamics of its coefficients. We conclude this section by deriving from it a deterministic equation for 
the EDF. The solution of this equation allows us to calculate the van Hoves function and the dynamic 
structure factor of the gas, which are obtained in Section 4. In Section 5 we show that when the long 
time limit of the EDF is evaluated as a function of the position, the contraction procedure and the 
non-Gaussian corrections introduced by solving the MSE, yield a behavior which reproduce the one 
proposed by NESM for the same model. We conclude with a discussion of the results of Section 5. 

II. SMOLUCHOWSKI EQUATION 

Following Ref. we consider a driven nonequilibrium system composed of regions (cells) of size / 
where fluctuations of 1/T may occur [iQ]. In this model system, it is assumed that the local temperature 
of a cell changes in a time scale, r, much larger than the relaxation time a region needs to reach local 
equilibrium. Within each cell a spherical test particle of radius a and mass unity, performs Brownian 
motion and then moves to another cell. Its velocity may be described by a linear Langevin equation of 
the form 

dv ^ , , , , , 

— = -jv + aL{t), (1) 

where 7 > is a friction constant, L{t) is Gaussian white noise and a describes the strength of the noise. 
For this model the inverse temperature /5 = is no longer constant, but fluctuates in space and 

time on the scales / and r, respectively. As a result, Brownian motion takes place on two time scales. 



one related to the dissipation of the kinetic energy through the Stokes friction coefficient 7 = Qirrja, 
where 1] is the shear viscosity of the host fluid, and the second one is associated with the presence 
of temperature fluctuations induced on the whole system by an external agent. These fluctuations 
originate in the environment of the Brownian particles (heat bath) and in some cases may therefore 
possess a complex dynamics. 

In local equilibrium, the probability of finding a Brownian particle in position r with instantaneous 
velocity u in the presence of the temperature fluctuations 6T, is described by the local equilibrium 
probability density Pie{r,u;6T). This quantity may be determined through the calculation of the 
amount of work exerted on a particle at temperature T in the presence of an external potential (pjr ) 
and of a local temperature fluctuation 6T{t) = T — Tb, where Te is the temperature of the bath 12|, 
[3, If the spontaneous thermal (equilibrium) fluctuations AT are such that AT <^ 6T, the 

expression for Pie is IGj 

Pi^{u) = Poe"^^ +Ve~~^^, (2) 

where m is the mass of the particle, = de/dTs is the specific heat at constant volume, e is the 
internal energy density and where the density fluctuations of the heat bath have been neglected . 

The dynamics of the Brownian particle may be described through the nonequilibrium probability 
density P(r, u, t; 5T). The time evolution equation for this quantity is derived from the entropy produc- 
tion dAs{t)/dt, where the nonequilibrium entropy increment As(t) is given by the generalized Gibbs 
entropy postulate 

As{t) = -ke I Pln^dfdu. (3) 

J Pie 

Since the probability conservation is expressed through the continuity equation 

?i^-l.(up)-i.m, (4) 



dt df du 

from Eqs. ([2]), ([3]) and (jl]) and by using the boundary conditions for the fluxes, we get an expression for 
the entropy production where fluxes and forces may be clearly identifled, [l^, [l^. By assuming 



linear relationships between fluxes and forces we flnally arrive at the Fokker-Planck equation 

(5) 



d p ^ d _ d 

dt df du 
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where we have taken the couphng coefficient for the external force d(j)/df as one, in accordance with 
the phenomenological description we are using [lol, [2^. It should be pointed out that in arriving at 
Eq. ([5j) we have neglected the contribution arising from the usual thermal fluctuations in Eq. ([2]), in 
accordance with the assumption AT <^ 6T. 

The Smoluchowski equation for the reduced probability distribution p(r, t) = J Pdu can be derived 
in the limit of long times, by calculating the evolution equations for the ffist three moments of P over 
M-space, namely, the mass density p(r, t) = J Pdu, the momentum density pv{f, t) and the stress tensor 
density. For sufficiently long times, one may assume that the stresses and momenta have relaxed, and 



by combining them we obtain a constitutive relation for the diffusion flow (see Ref. 16| for details) 



After substitution of this constitutive relation into the continuity equation we are lead to 

— = ■ fp— ^ + (6) 

dt dr \ df J dr^ ' 

where in the last term we have taken into account that T = Tb + ST and Eq. ([6]) contains the diffusion 
coefficient 

D = (kBTB/mP) (1 + 5T/Tb) , (7) 

which explicitly incorporates the effects of the externally induced fluctuations on the temperature of the 
bath. It is convenient to stress that equations similar to have been previously derived in the context 
of slow relaxation systems such as supercooled colloidal suspensions and granular systems, among others 



2l| . 221]. However, in these cases, the temperature T was assumed to be a decaying function of time, 
related with an activated process which controls the relaxation. In contrast, in the present work we 
considered the more general situation where the temperature fluctuation ST{t) is assumed to be a 
stochastic process, which converts Eq. ([6]) into a multiplicative stochastic partial differential equation 
(MSE). 

III. EFFECTIVE PROBABILITY DISTRIBUTION 

With the purpose of obtaining macroscopic local variables which eventually may be measured, such 
as the dynamic structure factor, the van Hove's function, or the local mass density, we contract the 
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description over velocities and ST{t). To this end we define the EDF as 

(p(f,t)) = [ P{f,u,t;6T)dud{6T). 



A time evolution equation for the EDF is obtained by deriving an approximate solution of Eq. ([6]) 
;hrough the use of a method proposed by van Kampen which yields an equation for the average (p(r, t)), 
isl . This approximation takes the form of a series expansion in powers of the Kubo number, ar^ 

which is assumed to be small. Here a is a parameter measuring the magnitude of the fluctuations in 

the coefficients of Eq. ([6]) and denotes their finite autocorrelation time [l5|. 

This derivation may be easily generalized to consider the case where the particle is also under the 

infiuence of a harmonic force derived form the potential (p = where c^o denotes its characteristic 

frequency. Let us assume that ST{t) is of the form 

6T{t)/TB = aCit), (9) 

with ({{) is an arbitrary stochastic process with zero mean and autocorrelation 

(C(t-r)C(t))=n(r). (10) 

Accordingly, we recast Eq. ([6]) as the following multiplicative stochastic equation for p{x,t), 

^ = V ■ (px) + rvV + aTCm'p. (11) 

where for convenience we have introduced the dimensionless variables x = a^^f, i = uj'^(3~^t and 
the operator V indicates differentiation with respect to x. The parameter F = kBTB/ma'^uj^ is the 
ratio between the thermal energy of the bath and the energy supplied by the external source, while 
maintaining the nonequilibrium state. By taking the Fourier transform over the space variables, Eq. 
(fTTll is rewritten as 

p{q,t) = -g ■ ||- r (g^ + aC(t)g') P, (12) 
where the caret stands for the Fourier transform of a quantity. 
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To implement van Kampen's method explicitly, it is convenient to identify the systematic and stochas- 
tic operators on the right hand side {r.h.s.) of Eq. (|T2|) as 

Ao = -q-^^-Tq\ (13) 

and 

Ai{t) = -aTC{t)q\ (14) 

Following van Kampen's method for times satisfying Tc ^ t ^ a^^, the average {p{q,t)) over the 
realizations of ((t) obeys by itself the approximate non-stochastic equation 

d{p) 



POO _ _ 

Ao + a^ rfr(Ai(t)e"^Mi(t-r)e-"^°) 
Jo 



(p). (15) 



dt 

After evaluating the action of the operators Aq and Ai, Eq. (ITSl) may be rewritten in the more compact 
form 

- ( q- (p) + rVC(t)g^(p), (16) 



dt 

where C{t) is defined in terms of the autocorrelation of the noise 



c{i)^ / dT{cmi-T)), (17) 







which is so far arbitrary. Note that the second term on the r.h.s. of Eq. (fT6|) incorporates in an average 
way the effects of the temperature fluctuations on the dynamics of the system. It also introduces time 
dependent non-Gaussian corrections and gives rise to non-stationary effects on the dynamics. 

If for simplicity in the discussion we consider the special case where the temperature fluctuations 
may be represented by a stationary stochastic process, C{t) = C(0) = Co, and the steady state solution 
of Eq. (fT6|) reduces to 

(p(g)) = Poe-^^%^«', (18) 

where po is a normalization factor. To comply with the approximation of van Kampen's method, this 
solution is expanded up to order and then the inverse Fourier transform of the resulting expression 
is taken. This leads to a solution which contains explicit non-Gaussian corrections to the steady state 
distribution function, 

(19) 



(p(x)) = -^e-^^" 



with f{x) = r — 6r ^x"^ + 3. For future convenience we expand the above expression in a Taylor 
series in x, 



(P(^)) 



po(27rr) 



-1/2 



* 1 I 

2r 



1 + ^a'Cofix) 



(20) 



Notice that a different (but similar) expansion can be carried out in terms of the well known formula 



e ^ ~ (^1 — , which leads 



to 



(p(x)),=.po(27rr)-i/Ml 



X 



2Tb. 



1 + ^a'Cofix) 



(21) 



IV. VAN HOVE'S FUNCTION AND STRUCTURE FACTOR 



Owing to its direct relation with experimental techniques, we shall now calculate the density corre- 



lation function 



231]. The space Fourier transform of this quantity is the van Hove function F, 



F{q,t) = {piq,t)p{qo,to)), 



(22) 



where % will be considered as a reference wave vector, to is an initial time and, as before, the bracket 
(.. ) denotes an average over the realizations of the external noise C(^)- 

The time evolution equation of this quantity is obtained by multiplying Eq. ( fT2l) by p(go,^o) and 
using Eqs. ( fT3lfT^ . yielding 



d 



dt 



-~F\q,t)= Ao + Ai(t) F\q,t) 



(23) 



where we have used the notation F*{q,t) = {p{q,t)p{qo,to)). Since Eq. ( |23!) is of the same type as Eq. 
(fT2|) . van Kampen's method also leads to the evolution equation 



(24) 



for F{q,t) defined by Eq. fl22|) . Again, for the special case where the temperature fluctuations may be 
represented by a stationary stochastic process, the steady state solution of Eq. (1241) has the form (|T8l) . 



FM = e-^" e^^'' . 



(25) 
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FIG. 1: a) F{q,t) vs. q. The red lines correspond to the case without temperature fluctuations. The remaining 
lines green to blue, are obtained from Eq. ([26|) for t = 0.1, 0.3, 0.6, 0.9, respectively, b) F{q,t) vs. t for fixed 
values of q = 0.5 (green) and q = 0.6 (blue). The value of the parameters are Co = 0.8, F = 1 and a = 0.8. 

Furthermore, using this result and the fact that Eq. ( 12^ is a linear hyperbolic partial differential 
equation which can be solved by using the method of the characteristics, we arrive at the following time 
dependent solution of Eq. fl2^ 

F(g, t)^(^l + j g-f (26) 

Again, we have only retained terms up to order to be consistent with the approximation used in van 
Kampen's method, and we have used the scaling variable z{q,t) = t — In Since the z dependence 
of the solution of Eq. (!24|) should decay with z, we have chosen the form e~^^ . Equation (!26l) actually 
imply that F{q,t) has a similar behavior as the probability density {p{q,t)). Figures la and lb show 
a plot of F{q, t) as a function of q and t for fixed values of t and q, respectively. They show that the 
presence of temperature fluctuations indeed modifles the behavior of F{q,t) with respect to the case 
where no intensive fluctuations are present. We have expressed the wave vector as g = qk, where k is 
a unit vector in the direction of q. It is important to point out that, in order to represent Eq. (126!) 
by Figs. 1 (and Eq. fl28|) through Figs. 2), we have chosen values for the coefficients F, a and Cq 
consistently with the approximations introduced in the previous section. This means that the physical 
parameters (T, u, (3) characterizing the system are determined. 

Fig. la shows that the external fluctuations induce signiflcant corrections to F{q,t) for values of 
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FIG. 2: a) S{q,u)) vs. g for = 1.8 and b) S{q,Lj) vs. q iov q = 1.6. The red lines correspond to the case 
without temperature fluctuations whereas the green line include temperature fluctuations for a = 0.8, Cq = 0.8 
and r = 1. 

q > 0.5 for various fixed time values (see caption). The red lines correspond to the case without 
(intensive) temperature fluctuations, whereas the green to blue lines correspond to the case when these 
fluctuations are present. This figure implies that intensive temperature fluctuations slow down the 
decay of the density- density correlation function. On the other hand. Fig. lb shows the corrections 
corresponding to the time dependence of van Hove's function for q = 0.5 (green) and q = 0.6 (blue). 
The Fourier's transform of F{q,t) with respect to time defines the dynamic structure factor S{q,Lj), 

S{q,uj) 



S{q,u) 



Juj ln\q\-^q'^{l-^a'^CW) 



Siqo,0) 

Expanding this function up to first order in its argument, we arrive at 



(27) 



S{q,Lu) = G ^{q,uj) 

where we have defined the inverse propagator 

F 



1 + §g^(l - ^a'Coq') 



(28) 



G-\q,u) 



1 + 2^'(1 - l^o^Coq") 



1 2 




+ 


oj In \q^\ 



(29) 



Eq. fl28l) explicitly shows that the presence of the intensive fluctuations indeed modify the behavior 
of the structure factor, as expected. The presence of powers of q higher than q^ suggests the existence 
of long-range correlations induced by intensive temperature fluctuations. In Figs. 2 we show the 
projections of the structure factor for fixed q (Fig. 2a) and for fixed uj (Fig. 2b). Note that in 
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Fig. 2b the green line is ~ 1.2-times larger than the red one for on = 0. If ~ 2.5 the green line 
has the same value as the red one and for ~ 10 the green line is ~ 0.6-times larger than the red 
one. Clearly, the presence of temperature fluctuations induces an anomalous behavior of the structure 
factor in the cj-fixed plane. It should be pointed out that in this case, our results are only valid for 

q < (ca2r)-i/2(i + Vl + W). 



V. COMPARISON WITH SUPERSTATISTICS 



In order to compare the results obtained in this work with those derived from the superstatistics 
approach reported in Refs. , we recall that in these references it is argued that the stochastic 

with a distributed /3, gives rise to the generalized canonical distributions of 



differential equation 
NESM Is I, 



Pnesm[x) 



1 I ^0 2 

1 H X 

n 



(30) 



where n is a positive integer and f3o = {(3)- Eq- ( |30i) contains two free parameters {I3q and n) that must 
be fitted in order to represent the function. Both parameters arise from the initial assumption that 
the inverse temperature is a distributed variable; however, the physical meaning of these parameters 
remains unknown. Moreover, it is interesting to notice that the form of Eq. fl30|) is independent of 
the nature of the variable of interest. For example, in Ref. 10|, a similar expression (with one more 
parameter) is used in order to fit the distribution of velocities in a fluid undergoing stationary turbulence. 



Figure 3 shows a comparison between our steady state distribution function, as given by Eq. (!20!) 
(blue line), and the above expression for Pnesm{x) (red line), for different values of the parameter n. 
The green line corresponds to Eq. (pTI) . The difference between these curves lies in the fact that the 
superstatistical method lOj, essentially analyzes the statistics of a system in local equilibrium, and 
considers the fluctuations of the intensive parameters to be of a static nature. On the contrary, Eq. 
( fT8|) represents the long time limit of the dynamic analysis that has taken into account the (irreversible) 
time dependent process associated with the intensive fluctuations. Moreover, the parameters entering 
into Eq. (120|) are determined by the nature of the system and must satisfy the conditions imposed by 
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FIG. 3: The blue line represents the stationary EDF as given by Eq. ()20p vs. for n = 18 and with a = 0.5, 
Co = 0.25, r = 0.13. The green line represents an expansion of the exponential in the stationary EDF by 
using Eq. ^ with 6 = 8 and a = 0.02, Cq = 0.12 and T = 0.12. 

the approximations of Sec. Ill, (tcQ; ^1). In order to find a good comparison between the red and 
blue lines, we have chosen the adequate order of the polynomial arising from the Taylor expansion of 
the exponential. 

It should be pointed out that in all cases the value of F controls the form of the peak, whereas the 
strength of the noise a and the noise correlation Cq control the tail of the distribution. The larger the 
values of a and Co, the larger the inflexion intermediate values of x. 



VI. DISCUSSION 



In this work we have analyzed the problem of Brownian motion introduced in Ref. [IQ] by using 
a mesoscopic approach {MNET). The time evolution of the nonequilibrium state of the system is 
described by means of a Fokker-Planck equation for the probability distribution in the phase space of 
the particles including local intensive parameter fluctuations. The basic difference with the point of 
view adopted in [lo| is that the dynamics of the (intensive) temperature fluctuations is explicitly taken 
into account by representing them as a stochastic process. 

An essential element in our approach is the contraction of the above description which yields a 
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stochastic multiplicative Smoluchowski equation for the (local) reduced probability density. The (ap- 
proximate) solution of this equation yields a deterministic equation for an effective probability distri- 
bution (EPD) valid for an arbitrary stochastic dynamics of the random coefficients. As a consequence, 
we calculated two measurable quantities, the van Hove's function and the dynamic structure factor of 
the system; however, we are not aware of any experimental results to compare with these predictions 
of our model. 

The most important conclusion that arises from our approach is that the long time limit of our 
EPD as a function of position, exhibits a similar behavior to the so called superstatistical distribution 
in nonextensive statistical mechanics {NESM), 1^. This result shows that the dynamics of intensive 
parameter fluctuations is fully consistent with Boltzmann-Gibbs statistics. The coarse-graining proce- 
dure and the time dependent non-Gaussian corrections introduced in solving the MSE, gave rise to 
non- stationary effects on the dynamics and generated those features of the EPD which are usually 
identified with the so called nonextensive effects of NESM. 

Our approach can be used for both, stationary and non-stationary fluctuations of the intensive 
parameters. However, it should be stressed that our approach and conclusions are only applicable to 
the model analyzed here and to the class of nonequilibrium states we have considered in the present 
work. If for other nonequilibrium systems with a complex dynamics similar conclusions can be drawn, 
is an open issue that remains to be assessed. 
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